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ly^ ■ Abstract Using relativistic tensor-bispinorial equations proposed in we solve 

O , the Kepler problem for a charged particle with arbitrary half-integer spin inter- 

^ ! acting with the Coulomb potential. 

X) '■ 
<D ■ 

ptn : 1 Introduction 

T— I ■ 

^ ■ Exactly solvable problems in quantum mechanics are quite important and il- 

;> , lustrative but rather rare (see, e.g., |2]). They can be described fully and in 

'nI" I a straightforward way free of various complications caused by the perturbation 

I method. The very existence of exact solutions of these problems is usually con- 

CN ■ nected with their non-trivial symmetries which are mostly of particular interest 

^ ■ by themselves In addition exact solutions form complete sets of functions which 

O ■ can be used to find solutions of other problems. 

r-| I One of the triumphs of the Dirac theory for the electron is that it solves the 

I I problem of electron motion in the Coulomb potential. The exact Sommerfeld 

formula for the related energy levels is one of the cornerstones of relativistic 
^ ■ quantum theory. However extension of these results to the case of charge particles 

>• . with higher spin appeared to be very complicated since even in the case of spin 

^ ! s = 1 the corresponding relativistic wave equation (namely, the Kemmer-Duffin 

^ I equation) is not free of inconsistences and predicts the orbital particle will fall 

■ - - ' down on the attracting center [3] , j3] . The other well known problems of RWE for 

particles with higher spins are: violation of causality ill-defined interaction 
with a constant and homogeneous external magnetic field jB], and so on (for 
details see P, [Zj). 

In recent paper P relativistic wave equations for particles of arbitrary half- 
integer spin were proposed which are free of inconsistencies typical for other wave 
equations for higher spins. They are causal and allow the correct value of the 
gyromagnetic ratio g = 2. In addition these equations have well defined quasi- 
relativistic limits which admit a good physical interpretation and describe the 
Pauh, spin-orbit and Darwin couplings. 

In the present paper we apply the tensor-spinorial equations pLj to solve the 
Coulomb problem for particles with arbitrary half-integer spin. 



2 Relativistic wave equations for a particle with 
arbitrary half-integer spin 

First, we briefly review the RWE for free particles with arbitrary half-integer spin 
s and mass m described in pOj. The corresponding wave function ^[/^i'^il[/^2;^2] - [/in!^n] 
is an irreducible tensor with respect to the complete Poincare group of rank 
2n = 2s — 1 antisymmetric w.r.t. permutations of indices in the square brackets 
and symmetric w.r.t. permutations of pairs of indices ^ hJ = 

1,2, ■■■ ,n. The irreducibility requirement means also that convolutions w.r.t. 
any pair of indices and cyclic permutations of any triplet of indices of the wave 
function reduce it to zero. In addition, components of ■?/;[A'i'^il[/^2i^2]---[Mni^n] g^j^g 

bispinors of rank 1. This means that the wave function has an additional spinorial 
index a (which we omit) running from 1 to 4. 

The equation of motion has the following form P] 

(1) 

Here 'f^ are Dirac matrices, = i^- and the symbol Tip denotes the sum over 
all possible permutations of subindices (2, ■ ■ ■ ,n) with 1. 

In addition, the wave function ■?^[A*i'^i][M2i^2]---[Atn'^n] i^as to satisfy the static con- 
straint PP 

7/.7i.V^f^''"''''''^'"'''"''"' = (2) 

which is necessary to reduce the number of independent components of the tensor- 
spinor from 16s to 4(2s + 1). In order to obtain a theoretically recognized 2(2s-|- 
l)-component wave function we impose on ijj[f^i'^i]if^2'^2]-[p.ni^n] additionally either 
Majorana condition or a parity-violating constraint (1 + i75)'?/;[/^i'^iH/^2!^2]---[/i7i!^n] _ g 

(or (1 - Z75)^[A'l'^l][/^2l'2]-[Mnl'n] ^ 0). 

Equations (^, (j2I) are manifestly invariant with respect to the complete 
Poincare group and admit the Lagrangian formulation. For the case of a charged 
particle interacting with an external electromagnetic field this equation is gener- 
alized to the following form pQ: 

(^'y^^n^ — m) -?^[A'i^i][A'2i'2]-[A'ni'„] 

_J_5]p (7A'l7''l — ^I'l^Ml) 7]-;^^^^[A(T][M2;^2][M3l'3]--[Mnl'n] 

(3) 

,2iekJ^ (l^y 'J _pao-^[Mli'i] [M2i^2]---[Mn!^n] ^ ' 
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Here ttx = px — eAx, Ax and F'^'^ are vector-potential and tensor of the external 
electromagnetic field, k is an arbitrary parameter and 



so that 

Equation (jHl) contains both the minimal and anomalous interaction of a par- 
ticle with an external field Moreover, putting k = 2s — 1 for the anomalous 
coupling constant k we get the "natural" value = 2 for the gyromagnetic ratio 

Solving equation © for ^[^i^iH/^^^^]-!^"^-] obtain 

^[^lll^l][^J.2l^2]■■■[^J.nIyn] _ J_^^^A^[Mll^l][At2l'2]---[MnI^n] 

and then by substituting this expression into Q we get the following second 
order equation: 



Here S^^iy are generators of the Lorentz group whose action on the antisymmetric 
tensor-spinor ^t;*ifi][A'2!^2]---[Mni^Ti] given by the following formula 

^flU^[tJ'l'^l][tJ-2l'2]--itJ.niyn] _ i_ j ^ [M1 ^^l] [M2 i^2] - ' [^n i^n] 

_gfJ.Ul^WtJ'l]W2tJ.2]---WntJ.n] _|_ ^i/i/i ^ ] [l^2M2] - ' Km™] " 



where 5^'^'^ is the metric tensor with signature (+, — , — , — ). 

In accordance with its definition tensor ^|^i^il['^2i^2] [f^nv„] ^j-q^j^^^qj-j^^ according 
the representation L)(s-l/2,0)®D(l/2,0)©L'(0,s-l/2)®L)(0, 1/2) = D{s, 0)® 
D{s — 1, 0) © -D(0, s) © -D(0, s — 1) of the Lorentz group. Moreover, condition (j21) 
reduces this representation to D{s, 0) © -D(0, s) whose generators without loss of 
generality can be expressed as 

Sab — ^abcSci S^a — i^Sai Ct, & = 1, 2, 3, 

where Sa are matrices forming a direct sum of two irreducible representation D{s) 
of algebra so(3), Eabc is totally antisymmetric unit tensor of rank 3 and e is an 
involutive matrix distinct from the unit one and commuting with Sa- This matrix 
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can be expressed via the Casimir operators Ci = S^^S^^ and C2 = \snfi\aS'^^ S^'^ 
of the Lorentz group: 

e = C,C,\ (7) 

Thus instead of (0) we can study the equivalent equation 

where \1/ is a 2{2s + l)-component spinor belonging to the space of irreducible 
representation D{s, 0) © -D(0, s) of the Lorentz group, Sa {a = 1,2, 3) are direct 
sums of two (2s + 1) x (2s + 1) matrices forming the irreducible representation 
D{s) of algebra so(3). The components of the related tensor ^^i^i][^2i'2]--[Mni^n] 
be expressed via components of spinor \1/ by using Clebsh-Gordon coefficients (see 
Section 7). 

Taking into account commutativity of matrix i with matrices Sa equation (jH)) 
can be decoupled to two subsystems 

where are eigenvectors of i corresponding to the eigenvalues e = ±1. 

Thus in spite of a rather complicated form of the first order equations (jH)) 
we received the second order equation (jSj) for ^|j^i''il['^2»^2] [fj-nfu] relation 
which expresses ^[A'i-i][M2-2]-[Mn<^„] ^^^^^ ^[Mi-i][M2i^2]-[Mn.„]_ Moreover, equa- 
tion (jSl) can be reduced to This equation has been already solved for particle 
interacting with a constant and homogeneous magnetic field We shall see 
that equations Q are very convenient and are easy to handle in the important 
case when the external field is generated by a point charge. 



3 Radial equations for the Coulomb problem 

Consider a charged particle with arbitrary half-integer spin s and electric charge 
e interacting with an external electromagnetic field. When this field is generated 
by a point charge Ze the related vector-potential has the form 

a 

A = 0, Ao = -, (10) 
r 

where r = xf + + x"^ and a = Ze^ . 

Now we shall use the reduced version of the equations of motion given by 
expressions (0) to Since both equations Q corresponding to e = 1 and to 
e = —\ lead to the same energy spectrum we shall consider only the case e = 1 
and omit index e of function in the formulae of this section which follow. 
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For the states with energy E the corresponding solutions ^' can be written as 



= exp{-iExo)ij{r), (11) 

where ipi^) is a (2s + l)-component function depending on spatial variables and 
satisfying the following second order equation 



- A + ika^^ ] tp. (12) 



Taking into account the rotational invariance of equation (fT2jl it is convenient 
to expand its solutions in terms of spherical spinors Qji,^: 

^ = Ur)^] ,n, (13) 

where fi^ ; ^ are orthonormalized joint eigenvectors of the following four commut- 
ing operators: of total angular momentum square J^, orbital momentum square 
L^, spin square 5*^ and of the third component of the total angular momentum 
J3, whose eigenvalues are j(j + 1), /(/ + 1), s(s + 1) and m respectively. Denoting 
I = j — X we receive 



m = —j, — j + 1, ■ ■ ■ , j, and A = —s, — s + 1, ■ ■ ■ — s + 2m 



SJ, 



where rrisj = s if s < j and rrisj = j if s > j. 

The expressions for spherical spinors via spherical functions are given in the 
Appendix. 

Thus we expand ip{'x) in accordance with formula (|T3|) where ^ are radial 
functions and summation over A is imposed which takes the values indicated in 
the above. We note that the action of the scalar matrix S • r to the spinors 
flj j-x m is well defined and given by the formula P 

S-r Qj j_x m = rKl{,Qj j_y m, (14) 

where K^-'^ are numerical coefficients whose values are presented in the Appendix. 

Substituting (|T3|l into (fT2|) . using (|Tl|l and the following representation for 
the Laplace operator A 

^r'l) - lA . (15) 



where is the square of the orbital momentum operator L = r x p, we receive 
the following equations for the radial functions 

F^A = (16) 
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where F is the second order differential operator 

^ = ^ + -^ + (^ + «/0 (17) 

and M is a matrix whose elements are 

Mav = A(A - 2j - l)5,v + igoiK'^,^,. (18) 

Formula (fTI)|) presents the equation for the radial wave function of a particle 
with arbitrary half-integer spin interacting with the Coulomb field. 



4 Energy spectrum 

Matrix M is normal^ i.e., it satisfies the condition MM^ = M'^M. Thus it is 
possible to diagonalize it using some invertible matrix V: 

M ^ M= UMU~\ Mxy = ^xx'^x (19) 

(where fx are eigenvalues of M) thus system (fT^ is reduced to the sequence of 
decoupled equations 



1 _ 

Fi^ = (no sum over A) (20) 



where is a A component of vector ^ = . 

The other interpretation of ^•^ is as follows, one. Let be an eigenvector of 
matrix M corresponding to the eigenvalue vx- Then it can be represented as 

^(■^) = l^uj^ (no sum over A), (21) 

where ujx is the eigenvector of M which does not depend on p and is a multiplier 
depending on p. 



Changing the variables r ^ p = 2\Jvn?' — E^r, ^ ^ f = y ^^23^^ equation 
pn|l is transformed to the well known form 



where 



13= , „ ■ kl^{2] + lf + iv,-ia'. (23) 

Notice that equation ()2Hj) (but with another meaning of parameters kx) appears 
in the non-relativistic Hydrogen system pLO, and can be easily integrated. For the 
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bound states, i.e., for > E"^, its solutions can be expressed via degenerated 
hypergeometric function JF(n, d, p) as 



/ = Cp"^ exp (-0 ^ + I, p) 



(24) 



where C is an integration constant. 

Solutions ()24|1 are bounded at infinity provided the argument n = ^^^^ — P is 
a non-positive integer, i.e., n = —n' = 0, —1, —2, ■ ■ ■ . Then from we obtain 
the possible values of energy for bound states: 



E = m\ l + I . (25) 

{{n' + 1/2 + kxf -a-'Y 

Here kx are parameters defined in expression ()2H|1 . where ux takes the values 
which coincide with the roots of the characteristic equation for matrix M: 

det(M - i^xl) = 0, (26) 

where / is the unit matrix of the appropriate dimension D = 2s + 1 for s < j or 
D = 2j + 1 for j < s. 

Thus we present the exact values of energy levels for the Coulomb system for 
the orbital particle having arbitrary half- integer spin. However, formulae ()25|) . 

include parameter ux defined with the help of the algebraic equation of 
order D which can be solved in radicals for j < 3/2 or s = 3/2. For other values 
of s and j the formula defines an algebraic equation whose order is larger 
than 4, which does not have exact analytic solutions. The related possible values 
of ux should be calculated numerically. 

In the following section we find analytic expressions for approximate solutions 
of (j26|l and expand energy levels (j25|l in th power series of (ga)"^- 



5 General discussion of energy spectrum 

First we note that if g = 2, and s = 1/2 then equation reduces to the exact 
Sommerfeld formula for energy levels of the Dirac particle in the Coulomb field. 
Indeed, using and solving ipHjl for s = 1/2, g = 2 we obtain 



1 
4 



2Ai 



A = ±- 



so that the related formula takes the following form 



E = m 



1 + 



a 



n' + 



a' 
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i.e., coincides with the Sommerfeld formula. This seems to be rather curious 
since the Dirac equation appears as a very particular case of equation Q in 
which wave function ^[m^^il-lMn^^n] j^g^g 

zero number of pairs of anticommuting 
indices [/^ji^j]. The second order equation (fT^ with g = 2 and S being matrices of 
spin 1/2 , appears in the Dirac theory too if we use potential (fTIH) and express the 
"small" components of the wave function in terms of "large" ones (this procedure 
is equivalent to our substitution Q). Thus the spectrum of energies for higher 
spin fermions (^3)) includes the spectrum of the Dirac particle as a particular 
case. 

In order to analyze spectrum of energies of the orbital fermion with arbitrary 
spin we shall look for approximate solutions of equation (pUj) . Considering a as 
a small parameter, using explicit expressions (fTHj) and (jSHj) for matrix M and 
applying the standard perturbation technique, we obtain analytic expressions 
for approximate solutions of this equation for arbitrary spin and total angular 
momenta: 

ux = X'- (2j + 1)A 

(agf f {X + s){2j-X~s + l){s-X + l){2j + s-X + 2) 
8 V (2j-2A + l)(2j-2A + 3)(j-A + l) (27) 

(g + A + l){s - X){2] - s - X){2] + g - A + 1) 
(j-A)(2j-2A + l)(2j-2A-l) 

Starting with and using (f?fj) we find approximate expressions for energy 
levels up to the terms of order a^: 

E _ o? 

m~ ~ 2^^8n^~ n^{2l + 1) 

g^c^^ ( {j -1 + s){j + 1 - s + l){l - J + s + l){j + 
M\2l+lY\ (/ + l)(2/ + 3) 

U-1 + S + l)(j + / - s){l-3 + s)U + 1 + S + 1)(1 - 

l{2l - 1) 

n = n' + / + 1 = 1, 2, ■ ■ ■ , / = 0, 1, ■ ■ ■ , n - 1. 

Formula (j28|l defines the fine structure of the energy spectrum for particle 
with arbitrary spin. We see that energy levels are labelled by three quantum 
numbers, i.e., by n, / and j. Any level of order corresponding to the main 
quantum number n is split to sublevels which have the order and describe 
the fine structure of the spectrum. The value of coincides with the number 
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l + s + 2) 

(28) 

)• 



of possible pairs i.e., 



for n < s 



n 



(s + i)(2n-s-i) f 



or n > s. 



We recall that for the Dirac particle the number of split sublevels is equal 
to n and any level with j 7^ n — 1/2 is double degenerated. This again is in 
accordance with our formula ()28p if we put s = 1/2 and g = 2 and find the 
related approximate ux from (j2Zj). 

The physical interpretation of spectrum is obtained by using the results 
of paper P where the Foldy-Wouthuysen reduction of equations was carried 
out. In addition to the rest energy term proportional to m the formula 

includes the non-relativistic Balmer term —-^ and the term ( -,,1 , -.s ) 



which represents the relativistic correction to the kinetic energy. The last terms 
in (j^H|) which are proportional to g'^a'^, represent the contribution caused by the 
spin-orbital, Darwin and quadrupole interactions. 

6 Energy levels for s = 3/2 

In the previous section we have defined the energy levels for a charged particle 
with arbitrary half-integer spin interacting with the Coulomb potential and briefly 
described the related wave functions. Here we consider in full detail the case 



First we notice that in the case s = 3/2 it is possible to solve the corresponding 
algebraic equation in radicals and to find the associated exact formulae ()25|) . 
We shall not present these rather cumbersome formulae here and restrict ourselves 
to analysis of approximate solutions (j28|l . 

The energy levels ()28|1 for s = 3/2 are reduced to the following form 



Here j is the total angular moment quantum number which can take the following 




s = 3/2. 




(29) 



values 




(30) 
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and Aj i are corrections to energy levels defined by j and /: 



5(2/ + 1)2 I (/ + l)(2/ + 3) 



(31) 



(j - ^ + i)(j + ^ - i)(^ - J + i)(j + I + |)(1 - Sm)\ 1 



/(2/-1) y 2/ + 1' 

Thus in accordance with (jHUjl . (jHH) we have 

3c/2(2/ + 5) . 
^'^^'/^ ' - 8(/ + l)(2/ + l)(2/ + 3) " 27TT ^""^ ^""^ ^' 

^2(2/2 _ 9/ _ 27) 1 



3c/2(2/ + 5) 


1 


8(/ + 1)(2/ + l)(2/ + 3) 


2/ + 1 


^2(2/2 - 9/ - 27) 


1 


8/(/ + 1)(2/ + l)(2/ + 3) 


2/ + 1 


^2(16 - 2/2 - 13/) 


1 


8/(/ + l)(2/ + l)(2/-l) 


2/ + 1 


3^2(3-2/) 1 


— for 


8/(2/ + 1)(2/ - 1) 2/ + 


1 



(32) 

A;-i/2 1 = ^, " ^,,;7^ ^, - -— for/>o, 



A/_3/2 / 



Formulae (j^ - ljH^ describe the energy spectrum of a particle of spin 3/2 
interacting with the field of a point charge up to order a^. We see that if / > 1 
then any level corresponding to fixed values of the quantum numbers n and / is 
split to four sublevels representing a fine structure of the spectrum. The value of 
the fine splitting decreases as the quantum numbers n and / increase. 

We notice that unlike the case s = 1/2 energy levels ()29|) are in general 
non-degenerate. The fine structure of spectrum (f^n|) - (p?^ however includes an 
arbitrary parameter g which can be interpreted as a gyromagnetic ratio of the 
system pQ. This parameter can be fixed by means of experimental or theoretical 
requirements and some degeneracy can appear for its special values. 

There are three " privileged" values of g namely the generally recognized value 
g = 2 (which is predicted by string theory and is in accordance with experimen- 
tal data [H]), g = 1/s (which naturally appears in relativistic models without 
anomalous interaction [12]) and g = \/2/s. For g ^ \f2j~s it is possible to show 
that the spectrum ()29|) - (jn^ is degeneracy free in general and for (7 = 2 or = 1/s 
even accidental degeneracy free. However for g = a/2/s this spectrum is doubly 
degenerate. Moreover, as we shall see, the origin of this degeneracy is similar to 
one observed in the Dirac-Coulomb system (cf. [T!^). 
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(33) 



Setting in (jH^ g = y | we obtain for s = 3/2: 

4/(/ + 2) + l 

A; ,3 7 = — -— -— for any 

'+2 ' 2(2/ + l)(/ + l)(2/ + 3) ^ ' 

4/^(3/ + 7) + 27(/ + l) 
'+^' 6/(/ + l)(2Z + l)(2/ + 3)' ' 

4/^(3/ + 2) + 7/ -16 
'-5' 6/(2/-l)(2/ + l)(/ + l)' ' 

2/(2/ + l)(2Z-l)' 

Comparing A;_,_3 ^ with A;_3 ^ in ()33|1 it is not difficult to observe a double 
degeneracy of the related energy levels, due to the fact that 

A,+ 3 , = A, 1 ,^1, 1^0. (34) 



Thus if the gyromagnetic ratio is equal to y | then like in the Dirac- Coulomb 
system the energy levels corresponding to maximal or minimal possible value of 
j for given / > are doubly degenerate. However, in some ways, this analogy is 
rather matter of convention since in the Dirac-Coulomb system all energy values 
(except the ground one) are doubly degenerate. In our system with spin 3/2 like 
in the Dirac-Coulomb one there are degenerate states with j = I ± s. But in 
addition for s = 3/2 there are states with j = / ± (s — 1) and with I = which 
are singlets. 

We note that analysis of energy spectrum ()28p presented in this section admits 
a straightforward extension to the case of any spin. In particular, our conclu- 
sion concerning degeneracy of the spectrum for special value g = -^/| of the 
gyromagnetic ratio is valid for any s. 



7 Explicit solutions of the relativistic wave equa- 
tion for s = I 

Following the procedure outlined in Sections 4 and 5 it is not difficult to find 
the explicit expressions for the wave function of the particle with spin 3/2 in the 
Coulomb field. 

We start with the basic equation 0. In the case s = 3/2 the related wave 
function is a tensor-bispinor ■ip^^^'^^ which has 24 components, and equation ^ is 
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reduced to the following form 

(35) 

(h«7.i^"^V^f' + i^.^V^t"' - i^.^f ') = 0, 

where (and in the following) 

^M=^[H^l^^^M-^^^[Aa]. (36) 

Solving RWE for tp^^'^^ we obtain equations ^ and for particle with 
spin 3/2: 

(^^^a _ ^2 _ i^^^F-) = 0, (37) 



= l^^^vrV!^"^!. (38) 
m 

Functions ijj^^'^^ and ^/'l^;'''^ both have 12 independent components and form 
carrier spaces for representations £'(3/2, 0) © D(0, 3/2) © D{l/2, 0) © D(0, 1/2) 
and -0(1, 1/2) © _D(l/2, 1) of the Lorentz group respectively. 

In addition, ip^^'^^ has to satisfy the condition (j21): 

-ff.lu^^^"^ = 0. (39) 

In accordance with its definition, tensor ip^^'^^ automatically satisfies the con- 
dition 'Jn'-fy'ip^!^'^^ = (even when condition ()39|1 for ip^i^'^^ is not imposed). Function 
ip\l!:'^^ satisfies condition '^^^'^^4'+'^^ = as well provided equation pTIjl is satisfied. 
Moreover, relation reduces the number of independent components of tp\[!^'^^ 
to 8 which form a carrier space of the representation D[3/2,0) © D(0, 3/2). 

The action of the generators of the Lorentz group on ip^^'^^ can be derived from 
equation (jHl) for n = 1 : 



In fact the same action is valid for any tensor-bispinor, e.g., for and ip^^'^^ 
For tp^^'^'^ the formula (pn|) can be specified taking into account equations (pUjl and 
(jSH). First we note that in accordance with (jHUI) 

V^f ] = l5°V75^f ] (a, 6, c= 1,2, 3) (41) 
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so that we can restrict ourselves to spatial components ■i/'l^*^^ of the tensor-spinor 
with a, 6 = 1, 2, 3 (we denote superindices running from 1 to 3 by Latin letters). 
Denoting 

and using the relations (jHBj) . (jH^ . PU]) and (PT|) we obtain the action of generators 
of the Lorentz group in the following form: 

gab^c ^ ^abkgk^c^ gOa^c ^ ij^S^T]'', (43) 

where 

Sa = Si^'^ + Si'\ Si^^ = -^Eabclblc (44) 

and S^^ are spin one matrices whose elements are {^S^^^ = iSabc- 

We note that like tensor-spinor vector-spinor 77" = r/^ has an implicit 

spinorial index a which we usually omit. Moreover, 7- matrices act on spinorial 
index of r]"- while matrices S^^ act on vector index a so that S^^ commute with 7^ 
by definition. Thus matrices Sa in (jl^ are sums of commuting spin 1/2 and spin 1 
matrices, and so that they can be reduced to the direct sums of spin 3/2 and spin 
1/2 matrices. Such reduction is easily handled my means of Wigner coefficients 
provided the corresponding third components of spin vectors are diagonal. 

If we choose the following realization of the Dirac matrices: 

70 = ( A ? ) ' 7a = ( ° a" ) ' ^75 = ( ^ ° ) (45) 



^ -1 y ' ' \aa y ' " V . 

where I and () are the 2x2 unit and zero matrices respectively and ctq are the Pauli 

(-) 

matrices, then both ^75 and 5*3^ are diagonal. Using the unitary transformation 
ly^b^b ^j^}^ transformation matrix 



(46) 




we diagonalize S^^ too since in this case 




Thus we pass from fl44|) to the following realization of spin matrices 

Sa - Sa = WSaW^ = ^(Ta + Si'\ (47) 
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Now we are ready to apply the Wigner coefficients C^^_^g^^^ to reduce (P7j) to 
a direct sum of spin 3/2 and spin 1/2 matrices. In our case Si = 1, S2 = 1/2, s = 
3/2,1/2 and the following notation is used: m = k, rrii = c, m2 = a. Then 
the basis \E'*'^ which corresponds to the completely reduced representation of spin 
matrices is connected with 77^ by means of the formula: 

^sk ^ Cll,_W,ar]a- (48) 

Condition Q suppresses the states with s = 1/2 so that \1'^'^ = 0. Then, 
denoting ^fa^ = v]/'^ and using ()48|1 we obtain: 

= C'f^a^-^S- (49) 

Thus relations ()4H1 . P2|l and P9|l explicitly show how functions ijj^^'^^ (which 
satisfy equations (jSZj)) are connected with solutions of the reduced equation 
(jHl). The inverse relations have the form 

In accordance with expressions solutions of the reduced equations 

(jH} corresponding to energy E = En>j ux can be written as 

*n'i ux = eXp{-iEn'j ^^Xo)^a6^| j^X (51) 

where uj\ are the eigenvectors of matrix M whose elements M\\i are given by 
formula fTHj) (see (jSZj)), and: 

6 = C7,(m2 - El^ ^J^r^e^'"'-''"'^ -'a^'^^ (-n', fc, + 1, 2^m^ - E^^ ^j) . 

(52) 

Substituting = \E'^,^- from (j5T|l into (j50|) we obtain the explicit expres- 
sions for (non-normalized) solutions of the equation 

To complete presentation of solutions of the reduced equation (jH)) for spin 

3 

s = 3/2 we present the explicit expressions of spherical spinors j_A m 
Appendix. 

8 Discussion 

In the paper the explicit solution of the relativistic quantum mechanical Kepler 
problem for orbital particle with arbitrary half-integer spin are presented. As a 
mathematical model of such a system we use the tensor-spinor relativistic wave 
equations (RWE) proposed in paper [T]. 
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Tensor-spinor formulation of RWE makes it possible to overcome the main fun- 
damental difficulties which appears in relativistic theory of particles with higher 
spins like causality violation and complex energies predicted for particles inter- 
acting with the constant magnetic field |P . In addition, this formulation is rather 
straightforward and simple which makes it possible to extend the known exact 
solutions for some special quantum mechanical problems described by the Dirac 
equation 2J to the case of an arbitrary half-integer spin. The problem of inter- 
action of arbitrary spin fermions with constant and homogeneous magnetic field 
was solved in paper iS]. 

In the present paper we solve a much more complicated problem: interaction 
of arbitrary spin fermion with the field of a point charge. In addition, tensor- 
spinor RWE admit exact solutions for the cases of interaction with the plane wave 
field, with constant crossed electric and magnetic field and some others. We plan 
to present solutions of the related problems in the future. 

In addition to the general treatment for an arbitrary spin we restrict in full 
details to case s = 3/2. However the results given in Sections 6 and 7 admit 
straightforward generalization to the case of s arbitrary. In particular, in analogy 
with ()41|) - (jM?|) the tensor-spinor ^||"i'^il[^2i'2]--[Mni^n] ^j^j-^ arbitrary n in view of 
equation (j2]) is effectively reduced to symmetric tensor fji^^t^2---iJ.n -^vith n indices 
running from 1 to 3, and then to a 2(2s -|- l)-component function ip satisfying 
equation (fT^ . 

We obtain the generalized Sommerfeld formula ()25|) for energy levels of par- 
ticle with arbitrary half-integer spin interacting with the Coulomb potential. In 
contrast with the formula generated by the Dirac equation our expression (j25|) 
includes parameter g whose value is not fixed a priori. In accordance with the 
analysis present in pp this parameter is associated with the gyromagnetic ratio 
of a particle described by equations (0), ©• 

Analyzing formula ()28|) for energy levels we conclude that in addition to the 
most popular values g = 1/s and g = 2 there exist one more intriguing value, 
namely g = \/2/s which corresponds to a specific degeneracy of the related 
energy spectrum. We notice that in the case s = 1/2 all mentioned privileged 
values of the gyromagnetic ratio coincide while for s > 1/2 the relation 



is valid. In other words the degeneracy related value of g lies between the recog- 
nized values 1/s and 2. 

One of us (AGN) is indebted to the Ministery of Education, Youth and Sports 
of the Czech republic for the support throw the grant number 1P2004LA211. 
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A Appendix 

Here we present some technical data used in the main text. 

Spherical spinors are eigenvectors of the commuting operators of squared total 
angular momentum J^, squared orbital momentum L^, squared spin and the 
third components J3 of total angular momentum, so that 



^^j I m 
^^j I m 



I ml 



Spherical spinors can be expressed via spherical functions Yim' 



Yi 



where C/™_^ ^ ^ are Wigner coefficients, 



^„ if) = ^exp(.rnv^)(-l)-^ i'^^mi^Y^ Pri-osie) 



Y, 



pl'"l(cos( 



{l+\m\y. 

— sini™i e- — ^(cos' e - ly, 

m (dcose)H^ 



Lf and 9 are polar and azimuthal angles of r. 

The result of action of matrix S ■ r to spherical spinors is given by relation 
(dH) where 0, [121 



= -1/2 (^A'A+l a^+A + «s+A+l) > 

^ / /i(2j + 1 - /i)(2g + 1 - /i)(2j + 2s - /i + 2) 
"'^ V (2s + 2j-2;u + l)(2s + 2j- 2/2 + 3) 



1/2 



(53) 



In particular, for s = 3/2, the spherical spinors can be represented as the 
following four-component rows: 



/ ^(j +m){j+m- l)(j + m - 2)Y^ \ 



3 



'3 3-1 rn 2y/j(j~l){2j~l) 



3 2 "" 2 

'3-1 rn- 5 

v/3(j2 - m2)(j - m - 1)F^._| „+i 



v/3(j2 - m2)(j + m - 1)F._8 



V V (i - m)(j - - l)(j -ra- 2)Yj_i „+3 / 



/ -V^U + m)(j + m - l)(j - m + !)!;._ 1 „_3 \ 
-(j + 1 - 3m)^ yj + m Yj_i „_i 
(j + 1 + 3m)Vj ^+1 



\ V3(j - m)(j - m - l)(j + m - l)rj._i „+| / 
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3 



j j+2 ™ 20(i+l)(2i+3) 



f \/3(j + m)(j - m + 2)(j - m + „_3 \ 

-(j + 3m) V j - m+T r^.^i „_i 
-(j-3m)v/j + m + iy^.^i 
V V3(j - m)(j + m + + m + 2)Y^^i ^+3 / 



3 

.^3 



/ -VU + - ^ + 2)(j - + 3)y^.+ 3 ^_3 
v/3(j + m + l)(j - m + l)(j - m + 2)F^.+ 3 „„i 
-v/3(j - m + + m + + m + 2)Y^^ 



V v^(j + m + + m + 2)(j + m + 3)Fj 
and the coefficients ()53|1 are reduced to the following form: 



^3 (j + l) 



^7 



( 2j-l)(2j + 3) 



, 2-' 



3j 



J + 1 



(54) 



The related matrix M whose elements Mw are given by formula p8|) takes 
the form: 



/ ¥ + 3j 



M 



iga _3j_ 

2 V i+1 







igo' I 3j 3 _|_ ■ 
2 V j+1 4 



iga / (2i-l)(2j+3) 
2 V 









iga / (2j-l)(2i+3) 
2 V 





J 



iga / 3(2+1) 



iga 1 3(j+l) 
2 V i 



3j 



(55) 



Eigenvalues v of matrix coincide with the roots of the characteristic 
equation for matrix which is an algebraic equation of order four: 



2zy4 _ i0i/3 + (f + b{gaf - 10j(j + 1)) 

+ (i8j(j + 1) - 1 - 'iigaf) y + 18j2(j + 1)2 

+ 1) - I (!)' - (9J(J + 1) - 48.75) (^7«)2 = 0. 



(56) 



The related eigenvectors ujx corresponding to eigenvalues z/a, A = 1, 2, 3, 4 can 
be observed in the following form 



/ 



(pa)V3K2.-l)(2,+3) ^^^^^^^^._3^ 



\ 



»(ga)V(2i-l)(2j+3)^ 



((15 + 12j-4z/a)(3 + 4j-4z/a) + 



15 + 12j-4z/a)(3-12j-4z/a) 



12(ga)2(j+l) 



(4z/A + 12j - 3) 



v 



iga 



V3(i+1) 



2VJ 



(15 + 12j-4z/a)(3 + 4j-4z/a) + 



12(ga)2(j+l) 



(57) 
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Relation (j37|) is used in formula (jHH)- 

Solving algebraic equation it is possible to find exact values of ux. We 
shall not present the related cumbersome formulae here whose expansion in power 
series of {gaY is given in formula ()27|1 . 
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